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Summ£iry. The paper presents results for deriving closed-form analytic so- 
lutions of the non-relativistic linear perturbation equations, which govern the 
evolution of inhomogeneitics in a homogeneous spatially flat multicomponent 
cosmological model. Mathematical methods to derive computable forms of the 
perturbations are outlined. 
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1. Introduction 



This paper deals with solutions of a system of differential equations describ- 
ing the evolution of density perturbations in the Universe. Before being more 
specific about the differential equations, wc look at previous work setting the 
physical environment for deriving such equations. 

To formulate a quantitative picture of the evolution of perturbations outside 
the horizon it is necessary to solve the perturbed EINSTEIN equations. The 
density contrast d{x) is not a gauge invariant quantity. For sub-horizon-sized 
perturbations, aX = Xphys ^ H~^, this fact is of little consequence, as a NEW- 
TONian analysis is sufficient; a denotes the scale factor and H the HUBBLE 
parameter, respectively. For super-horizon-sized perturbations, Xphys ^ H~^, 
the gauge transformation for S{x) must be taken into account. In this case a 
full general relativistic treatment is required. So far two approaches are avail- 
able for finally deriving the differential equation governing the growth or decay 
of gravitational instabilities in an expanding Universe. The metric perturba- 
tion approach was invented by LIFSHITZ (1946) with his derivation of the 
synchronous equations. This approach was later used by BARDEEN (1980) 
to derive the comoving equations. HAWKING (1966) initiated an alternative 
approach to follow the evolution of the density perturbations in the expanding 
Universe by employing the general relativistic fluid flow equations. This method 
was used by OLSON (1976) to derive equations which are equivalent to the syn- 
chronous equations, and the same method was further developed by LYTH and 
MUKHERJEE (1988) to provide the comoving equations. 

In order to split up the energy density at a given point in spacetime into an 
average value plus a perturbation, one must specify a spacelikc hypcrsurface on 
which the averaging is to be performed. The choice of coordinates in order to 
define the energy perturbation throughout spacetime is referred to as a choice 
of gauge (BARDEEN 1980). BARDEEN proposed to use the comoving gauge 
to deal with the general relativistic equations. Another choice is to use the 
synchronous gauge (PEEBLES 1980) which, however, is besot by the problem 
of introducing arbitrary "extra gauge mode" solutions. Actually there is an 
infinity of synchronous gauges. 

A choice of LIFSHITZ's or HAWKING's approach and synchronous or co- 
moving gauge having been made sets the framework for deriving differential 
equations describing the time evolution of the energy density perturbation. Here 
again one has the choice to subscribe to the perfect fluid approximation or the 
kinetic description (PEEBLES 1980). In the perfect fluid approximation, which 
is more physically transparent, the time evolution of the energy density pertur- 
bation is specified by a differential equation involving the pressure perturbation 
as an additional unknown. The equation is of second order in the comoving 
gauge, assuming that the pressure perturbation is known, and it has two so- 
lutions revealing a growing and decaying mode (PEEBLES 1980). Since the 
equation governing the amplitude S{t) is of second order, there are two solu- 
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tions. A given perturbation is expressed as a linear combination of d+{t) and 
S-{t). At late times, only the projection onto the growing mode may be impor- 
tant. Physically, the decaying mode corresponds to a perturbation with initial 
overdensity and velocity arranged so that the initial velocity perturbation even- 
tually "undoes" the density perturbation. In synchronous gauge the equation 
is of third or higher order corresponding to the extra freedom of choosing the 
respective gauge. 

Many efforts have been made to derive approximate solutions of the equa- 
tions of relativistic and non-relativistic perturbation theory, usually valid for 
early time large-scale and for late time small- scale inhomogeneities, respec- 
tively, the primary purpose of this paper is to derive closed-form analytic so- 
lutions of the non- relativistic linear perturbation equations, which govern the 
evolution of inhomogeneities in homogeneous spatially flat multicomponent cos- 
mological models. The mathematical method employed here can also be used 
to tackle the relativistic version of the equations involved, which will be done in 
a subsequent paper. These closed-form solutions are valid for irregularities on 
scales smaller than the horizon. They may be used for analytical interpolution 
between known expressions for short and long wavelength perturbations. In 
Section 2 we present the physical parameters relevant to the problem and the 
fundamental differential equations for arbitrary polytropic index 7j and a radia- 
tion or matter dominated Universe, respectively. Solutions of this fundamental 
equation are provided in all subsequent Sections 3,4, and 5 for relevant param- 
eter sets including the polytropic index 7^ and the expansion law index ry. To 
catalogue these solutions in closed-form, in all Secitons the theory of MEIJER's 
G- function will be employed. The following results are a continuation of the 
ones by HAUBOLD, MATHAI, and MUECKET (1991) MATHAI (1989), and 
NURGALIEV (1986). 

2. Some Parameter of Physical Significance 

The growth of density inhomogeneities can begin as soon as the Universe 
is matter-dominated. This would be also the case in a Universe dominated 
by non-interacting relic WIMPs. Baryonic inhomogeneities cannot begin to 
grow until after decoupling because until then, baryons are tightly coupled to 
the photons. After decoupling, when the Universe is matter-dominated and 
baryons are free of the pressure support provided by photons, density inhomo- 
geneities in the baryons and any other matter components can grow. Actually 
the time of matter-radiation equality is the initial epoch for structure forma- 
tion. In order to fill in the details of structure formation one needs "initial 
data" for that epoch. The initial data required include, among others, (1) the 
total amount of non-rclativistic matter in the Universe, quantified by Qq, and 
(2) the composition of the Universe, as quantified by the fraction of critical 
density, fij = pi/pc, contributed by various components of primordial density 
perturbations (i= baryons, WIMPs, relativistic particles, etc.). Here the critical 
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density pc is the total matter density of the EINSTEIN-DE SITTER Universe. 
Speculations about the earliest history of the Universe have provided hints as 
to the appropriate initial data: fio = 1 from inflation; 0.015 < D.b < 0.15 
and ^wiMP ~ 0.9 from inflation, primordial nucleosynthesis, and dynamical 
arguments. In the following we assume that ^ Oj = 1, = const.. The cosmo- 
logical medium is considered to be unbounded and that there may be a uniform 
background of relativistic matter. It is further assumed that the matter is only 
slightly perturbed from the background cosmological model. This assumption 
may give a good description of the behavior of matter on large scales even when 
there may be strongly nonlinear clustering on small scales (primordial objects). 
Also assumed is that matter can be approximated as an ideal fluid with pressure 
a function of density alone. It consists of i components having the densities 
and velocities of sound 01 — dPi/dpi oc p'^^^^ , when an equation of state 

Pi cc p^- (1) 

has been taken into account. The i components of the medium are interrelated 
through NEWTON's field equation A(j) = ATiGj^Pi^ containing the combined 
density Yl Pi of all components. Superimposed upon an initially homogeneous 
and stationary mass distribution shall be a small perturbation, represented by 
a sum of plane waves 

5 = ^ = 5^(t)e»fc-, (2) 
P 

where A = 2'Ka/k defines the wave number k and A is the proper wavelength. In 
the linear approach the system of second order differential equations describing 
the evolution of the perturbation in the non-relativistic component i is 

where for all perturbations the same wave number k is used. This equation is 

valid for all sub-horizon-sized perturbations in any non-rclativistic species, so 
long as the usual FRIEDMANN equation for the expansion rate is used. Al- 
though the following investigation of closed-form solutions of the fundamental 
differential equation governing the evolution of inhomogcncitics in a multicom- 
ponent cosmological medium is quite general, it will be presented within the 
context of the inflationary scenario. Therefore in the following we consider an 
EINSTEIN-DE SITTER Universe with zero cosmological constant. After infla- 
tion the FRIEDMANN equation for the cosmological evolution reduces to 

(^)^ = ^f^ = ^p. (4) 

The continuity equation gives p oc and from eq. (4) one has a oc t^/^ 
and Pi = Q.il&'KGt^. During the expansion the HUBBLE parameter changes as 
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H = rjt ^, where the expansion law index is r/ = ^ in the radiation-dominated 
epoch and 77 = | in the matter-dominated epoch respectively. The wave number 
is proportional to so that k"^ 0f = kft'^'-^~'^~'*^^ in cq. (3), where now the 
constants ki come from both the wave vector and the velocity of sound. If fcj = 
the sound velocity 0i equals zero, in which case the adiabatic index 7j loses 
its sense. Defining the parameter 04 = 2(2 — 1] — 7;) that absorbs the adiabatic 
index as well as the type of expansion law we can write for eq. (3) using eq. 
(4): 
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t^Ut) + 2r]tS,it) + kfe^d, = ^y2 ^jSjit), i = 1, . . . , m, (5) 

and dots denote derivatives with respect to time. Wc introduce a time oper- 
ator A = td/dt and change the dependent variable, the density perturbation 
Si{t), dealing with the function $i instead of Si by setting Si{t) = The 
equation (5) for the function $i is then given by 

2 



where 



= t-'^Si, bi = kfe^ - a\ ai = 2{2-r]- 7^, a = -( ) . (6) 



Observing that Sfij = 1 and operating on both sides of the differential equation 
for the $i by we obtain the fundamental equation 

2 2 " 

A''<i>i + A\bi^i)--{A^^i + bi^i) = --J2bj^j'^j, i=l,...m, (7) 

(Mathai, 1989). As indicated above for r] the values | and ^ are significant and 
in what follows it will be assumed that 2 > 7^ > |. We have chosen the range of 
values of 7^ for physical as well as mathematical reasons as will be evident later 
in the analysis of eq. (7). For some of these parameter values we will consider 
a multicomponent medium. Consider the special case 

61 = fc^-^-fcil^, (8) 

62 = bs = ---=bm = b = kH"-^^^l^ 

4 

of the fundamental equation (7). This is the case considered in Haubold, Mathai 
and Muecket (1991) in detail. In the present discussion we will use the same 
notations as in Haubold, Mathai and Muecket (1991). Let fc = in (8). This is 
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case 4.1 of Haubold, Mathai and Muecket (1991). For this case the parameters 
are the following: 



K,b*^ = ± 



(2r/-l ) 
4af 



2^2 



,bl,bl = ± 



1 



(2^-^ 2 
4 3 



4 3 3 ^ 



(9) 



Table 2.1 gives these parameters for 77 = |, i and 7j = |,l,|,|,2. 
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When evaluating the solutions explicitly one needs the parameter differences. 
These are given in Table 2.2. 
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3. Solution when ai = 

Note that when ^7 = | and 7i = | one can have ai = in equation (6). The 
exact solution given in equation (5.3) of Haubold, Mathai and Muecket (1991) 

is for the situation ^ 0. Hence this case needs separate discussion. When 
ai = the fundamental equation (7) reduces to the form 



(2r?-l)2 



(2r;-l)4 2(277-1)2 



+ '.3"^ 4 3 



16 



$1 = 0. 



(10) 



In this case the general solution is of the form 



(11) 



where ci, C2, C3, C4 are arbitrary constants and di,d2,d3, are the roots of the 
equation 



(2r? - 1) 



kl 



+ 



2 3 

(277-1)^ 2 (2r?-l)2 
3 



16 3 4 
They can be seen to be the following: 



(12) 



2 (2^ 2 

3 4 3 



kl 



0. 



(2r/-l)^ ^ 1 kl 



(13) 



dz,di = ± 



(2r/-l)2 



1 fc? 



1 fc? 



(14) 



4. Solution when the parameters differ by 
integers 

The general solution given in Haubold, Mathai and Muecket (1991), equa- 
tions (5.2) and (5.3), arc for finite values of t and for the cases that the 6*'s 
do not differ by integers. From Table 2.2 note that at two points, namely 
[rj = |,7i = |) and (77 = |,7i = 1) one has 6* - 63 = -1 and 6| - 62 = -1- 
This means that Gi and of (5.3) of Haubold, Mathai and Muecket (1991) 
need modifications. At all other points the Gj's remain the same. 
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4.1. Modifications of the G',s for the cases (ry = i,1i = 1) 



3' 

3 ' '* 3'^ 
Consider, for x = kft"^ /af^ai ^ 0, 



and (i] = |,7, - 



= G^^{x^^-l;}) (15) 
1 /■ r(i + s)r{-al - s)r(-a; - s)x-^ds 



^^^JL r(|-s)r(-i-s)r(|-s) 

Note that a zero coming from r(— ^ — s) at s = — ^ coincides with a pole coming 
from r(i + s) at s = -i. This can be removed by rewriting as follows: 

1 f r(| + s)r{-at - s)r(-aS - s)x-^ds 



2^^JL r(|-s)r(|-s)r(|-s) 

Evaluating (16) as the sum of the residues at the poles of r(| + s) one has 



,r( -af + f)r(-a^ + |) 

r(f)r(2)r(M) 



^ / * 5 , 5 10 ^ 14 
X 2i^3(-at + ^, -^2 + 4; — , 2, — ; -x). 

Evaluating G4 also the same way one has the following: 

G, = Gj;^ i:[t3'|+^) (18) 

1 /• r(-| + s)T{-al - s)T{-al - s) _ 

-X ds 



^ttiJl r(|-5)r(f-s)r(-i-s) 



'^^Ul r(|-,s)r(|-.)r(-i-5) 

_ -i n-i-ai)n-\-a*2) 

r(i)r(2)r(-i) 

X 2F3i---al,--~a*2;-,2,--;-x). 

Thus the complete solution in this case is of the form 

$1 = ciGi + C2G2 + C3G3 + C4G4, (19) 

where ci, C2, C3, C4 arc arbitrary constants, Gi and G4 arc given in (17) and (18) 
respectively and G2 and G3 are given in equation (5.3) of Haubold, Mathai and 
Muecket (1991). 
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5. Solution near oo for the parameter values of 
table 2.1 

Except for the points {r^^^i) = (|, |), (|, 2) in all other cases some of the 
parameters differ by integers. In five cases of Table 2.2 the poles in the integrand 
can be up to order 2 and in two cases the poles can be of order up to 3. The 
general solution near oo is of the form 

$1 = fiFi + /2F2 + /3F3 + /4F4, (20) 

where /i, /2, /s, fi are arbitrary constants and Fj are given in the equation (5.6) 
of Haubold, Mathai and Muecket (1991). For example, 

Fi = Gti {x , X = kle^/al, ai ^ 0. (21) 

F2 is available from Fi by interchanging al and F3 and F4 are available 
from Fi by replacing x by .xe*^ and .Te~*'^ respectively. We will evaluate Fi for 
the parameter values which differ by integers. 

5.1. Fi for (,y=l,7, = |) 
In this case one has 

F =J- f r(5)r(5)r(V6 + 5)r(-V6 + s)r(-at - s) 

^ 2m J L T{l + al + s) ■ K ) 

Note that ats = 0,— 1,— 2,... the integrand has poles of order 2 each. All other 
poles are of order 1 each. Thus 

Fi=Hi+H2 + H3, (23) 

where H2 and H3 are the sums of the residues at the poles of r(v^ + s) and 
r{—^/Q + s) respectively and Hi is the sum of the residues at the poles of r^(s). 
For all the cases of simple poles including H2 and Hs one has the result as 
follows: 



n-al + b*) 



J I k-i \ k jj\ r(l + -b*) ^ ^ 

X 2Fs{-al + b*,-a*2 + b*;l - b^ + b*, ,1 - bl + b*; -x), 

where 11' denotes the absence of the gamma r(&^ — 6^) and # denotes the 
absence of the element 1 — bj — bj. Thus one has, 

H2 = .v^r^(-76)r(-2V6) ^|-"^ + ^ (25) 

r(l + a2 - V6) 

X 2F3{-al + \/6, -a^ + \/6; 1 + \/6, 1 + \/6, 1 + 2\/6; -x) 
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and 

H, = .-v^r^(x/6)r(2x/6) (26) 

r(l + a^ + V6) 

X 2F3{-al - V6,-a*2 - V6;l - V6,l - V6,l - 2^6; -x). 

Now consider the evaluation of Hi. Note that for all parameter combinations 
(??,7i) = we have bl = = b^. Hence we can 

write Hi in general terms as follows, 

Hi = ±f rHs)m+fm+m-at-s)^_.^^_ 

2m J L T{l+al + s) ^ ' 

The general technique of evaluating such integrals is available in Mathai and 
Saxena (1973). The solution is the following 



Hi 



v=0 



^ (Inx) {2i}{v + 1) + V(6S -y)+ ^{bl - v) (28) 



+ij{-al^v)~ij(l + al-v)}] 



where the 63 and 6| are given in Table 2.1 for the specific parameter combina- 
tions and tp{) is a psi function. Note that the first part in (27) can be written 
as a 2-F3(— a*, — 1; 1 ~ ^31 1 " ^4! ~^)- 
5.2. Fi for {rj = |, 7i = 1) and (?^ = |, 7^ = |) 
In these two cases Fi is of the following form: 

p ^ J_ /• r(i + 5)r(-i + g)r(| + g)r(-| + 5)r(-at - s)x-^ds 

' 27rzA r(l + a* + s) ^ ^ 

1 /• r2(| + .)r2(-i + .)r(-.i-.)^_^^^ 



STTiA (s + i)(s-f)r(l + a^ + s) 

There are poles of order one each at s = —j and s = f respectively and poles 

of order 2 each at s = — | — z/, s = — A, = 0, 1, . . . , A = 0, 1, . . . respectively. 
Thus Fi is evalueated as the sum of the residues at all these poles. Let us denote 
these residues by Ri,R2,R3 and R4 respectively. Then 

Fi =i?i+i?2 + i?3 + i?4, (30) 

where 

Ri = residue at s = — - (31) 

4 
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^ r^(-i)r(-ai + 1) , _ 2 rH-m-< + 1) i. 
(i)r(| + a^) 3 r(f + a^) 

5 

R2 = residue at s = - 

_ 2 r^(|)r(-at + |) 

3 r(| + a^) 

°° r 3 5 

^3 = a;4^<^-lnx + 2V'(i^ + l) + 2V'(---i')+V(-aI + ^+i^) 
-V- 1 + a2 - - - + — - + jr—-, \ 

r2(-|-t.)r(~at + | + z.) X- 

^r(l + a^~|-^.)(l + j.)(| + ^.) 

1 °° r 3 

^(i_^)(_l_^)r(| + a^-j.) 

As before Fi is available from the F\ of (30) by interchanging a\ and F3 
and F4 are available from Fi replacing x by e'^^x and e~''^x respectively. This 
completes the evaluation of $1 for all finite values as well as for values near 00 
for all the parameter combinations given in Tables 2.1 and 2.2. 

6. Conclusion 

We have presented closed-form solutions of the non-relativistic linear pertur- 
bation equations which govern the evolution of inhomogeneities in a spatially 
flat multicomponent cosmological medium. The general solutions are catalogued 
according to the polytropic index 7j and the expansion law index 77 of the multi- 
component medium . All general solutions are expressed in terms of MEIJER's 
G-function and their MELLIN-BARNES integral representation. The proper 
use of this function simplifies the derivation of solutions of the linear pertur- 
bation equations and opens ways for its numerical computation. In this regard 
the paper leaves room for further work. 
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